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Abstract. A scmilincar reaction-diffusion two-point boundary value problem, whose second- 
order derivative is multiplied by a small positive parameter e 2 , is considered. It can have multiple 
solutions. An asymptotic expansion is constructed for a solution that has an interior layer. Further 
properties are then established for a perturbation of this expansion. These are used in [6] to obtain 
discrete sub-solutions and super-solutions for certain finite difference methods described there, and 
in this way yield convergence results for those methods. 
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1. Introduction. We are interested in interior-layer solutions of the singularly 
perturbed semilinear reaction-diffusion boundary- value problem 

(1.1a) Fu(x) = -e 2 u"(x) + b(x,u) = for x£ (0,1), 

(1.1b) u(0)=g , u(l)= 9l , 

where e is a small positive parameter, b G C°°([0, 1] x R), and go and g\ are given 
constants. 

Under the hypotheses that are stated below, this problem can have multiple solu- 
tions that exhibit interior-layer behaviour. A companion paper [S] discusses numerical 
methods for its solution. In this present report we shall present some of the details 
that are omitted from [6]. 

The reduced problem of (1.1) is defined by formally setting e = in (|l.la[) . viz., 

(1.2) b{x,ip) = for x e (0,1). 

Assume that the reduced problem (|1.2p has three simple roots ip = <fk e C°°[0, 1] for 
k = 0,1,2: 

(Al) b(x,(p k (x)) =0 for k = 0,1,2 and x 6 [0,1] 

where 



tpi{x) < ipo(x) < (f2(x) for x e [0, 1] 

and there is no other solution of (|1.2|l between (fi and if 2 ■ 



(A2) 



Here and subsequently, numbering such as (Al) indicates an assumption that holds 
true throughout the paper. Assume also that 

(A3) b u (x, <p k (x)) > 7 2 > for k = 1, 2 and x £ [0, 1] 
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but 

(AA) b u (x,ip ) < for x G [0, 1]. 

Assumption (A3) says that <pi(x) and <f2(x) are stable reduced solutions, i.e., one 
may have a solution u of that is very close to either </ji or <^2 on some subdomain 
of (0, 1). Assumption (AA) means that the solution <po(x) is unstable: no solution of 
(jl.ip lies close to ipo on any subdomain of (0,1). Under the hypotheses (Al)-(AA), 
the equation (|l.la|) is often described as bistable. 

Our further assumption is that the equation f^ffi b(x, v) dv = has a solution 

x = to such that [ f£*(*\ b(x, v) dv\ \ x _ t ^ 0, i.e., this root is simple. As in many 
asymptotic analysis papers, we also assume that the value of this derivative is negative, 
since this sign corresponds to the Lyapunov stability of an interior-layer solution u(x) 
of that switches from ipi to if2 when u is regarded as a steady-state solution 
of the time-dependent parabolic problem vt — e v xx + b(x, v) = (see [1] Section 7, 
Remark 3]; if instead the derivative were positive, this would correspond to Lyapunov 
stability of an interior-layer solution that switches from (f2 to ipi). By Assumption 
(Al) these hypotheses on the integrals of b are equivalent to the assumptions 

(Ab) / b(t o ,v)dv = and / b x (t , v) dv = -Cj < 0. 

J<pi(to) J<pi(f ) 

Similar conditions are assumed in [5J §4.15.4], [SJ §2.3.2] and also in [5J [7] for an 
analogous two-dimensional problem and [3] for a analogous system of equations. 

Remark 1.1. Assumption (A2) can be relaxed to allow other roots of \1.2\i 
between ipi and ip2 provided that J v ^^ tg j>(to, s) ds > for all v G (tpi(to), <P2(to)). 
Note that this inequality follows immediately from (A1)-(A5) if ipo is the only reduced 
solution between ipi and (f2- 

The solutions tpi and y>2 of p.2p do not in general satisfy cither of the boundary 
conditions in (|l.lb|) . In order to focus on interior layers, we exclude boundary layers 
by assuming that 

(A6) Vi(0) = 90, <f*Q-) = 9i, Vi(0) = Va(l) = 0. 

Under Assumptions (Al)-(A&), the problem (|1.1|) has a solution that, roughly 
speaking, lies in the neighbourhood of f\(x) and (p2(x) for x G [0,to) and x G (to, 1] 
respectively (see [6j Corollary 6.7]). Near x = to the solution switches from ipi to if2, 
which results in an interior transition layer of width 0(e \ lne|). 

The structure of the paper is as follows. In Section [2] an asymptotic expansion 
of a interior-layer solution of (jl.ip is constructed; this analysis draws on ideas of 
O El 13 E] • A modified version of this asymptotic expansion, related to [7], is examined 
in Section [3] This modified expansion is used in [6] to construct discrete sub-solutions 
and super-solutions for the numerical methods used in that paper to solve (jl.ll) . 

Notation. Throughout the paper, C,C',C and C' , sometimes subscripted, denote 
generic positive constants that are independent of e and of the mesh; furthermore, C 
and C are taken sufficiently large where this property is needed. These constants may 
take different values in different places. Notation such as / = 0(z) means |/| < Cz 
for some C . 
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2. Asymptotic expansion for the continuous problem. The point i £ 
(0, 1) is fixed by Assumption (Ab). Define the stretched variable £ by 

£ := (x- t )/e. 

Then a standard calculation shows that the zero-order interior-layer term Vo(£) of the 
asymptotic expansion of u is given by a solution of the following problem: 

(2.1a) - -^V o + b{t o ,V o ) = for £ e R, V {-<x>) = <pi{t ), Vb(oo) = ^(t ). 

We shall see shortly that (|2.1a[) has a solution Vb(£), but this solution is not unique 
as Vo(£ ± C) is also a solution for any constant C. Once we know that Vq exists and 
is a strictly increasing function, consider a specific solution Vq of (|2.1a[) subject to the 
parametrization 

(2.1b) V (0) = Mto). 

One might expect u(x) = tpo(to) to hold at x = to and thus the interior layer to be 
described by Vb(£). It is not the case, however; as we shall see below, u(x) = ip(to) at 
x = t = t,Q + eti+e 2 t2 + ■ ■ ■ , and the interior layer is described by Vo(£ — t\ — et 2 — •••). 
Here t±, t 2 , . . . are independent of e and can be found when constructing an asymptotic 
expansion of u, in particular, the values of t\ and t 2 are specified in the proof of 
Lemma 12.41 In our analysis, we take t = t + et\ + s 2 t2, skipping higher-order terms, 
and invoke a perturbed version of Vo(£ — ti — et^) defined by 

(2.1c) V (Z;p) = V (Z-t 1 +p), i 1 = h+Et 2 . 

Here the parameter p satisfies \p\ < p* for any fixed positive constant p* , but will 
typically take very small values. 

Lemma 2.1. Set j 2 :— mim^i^ b u (to, ipk(to)) > 7 2 ; where 7 > is from (A3). 
For any constant t\ and all \p\ < p* , there exist unique monotone solutions Vb(£) and 
Vq(^p) of WJa\) that satisfy (KTc\) , Huty . Furthermore, V and V are in C°°{R), 

(2.2) X(0 ■= iMO > 0, Xfop) ■= iVo&p) > for^ER. 
For any arbitrarily small but fixed X € (0,7), there is a constant C\ such that 

(2.3) m + xfop) < Cxe-^W for £ G R, \p\ < p*. 
There are constants C and C" such that for all \p\ < p* one has 

(2.4) C'x<V -Mto) <C" X for£<0, C' X < p 2 (to) ~V < C" X for £ > 0. 

Proof. In view of (A1)-(A5), these properties follow from the proof of [H Lemma 
2.1] or a slight extension of the proof of Lemma 2.1] using phase-plane analysis. □ 

An interior-layer solution u of problem (|1.1[) can be regarded as having a boundary 
layer at to on each of the subintervals [0,io] an d [to, 1]. Therefore we shall construct 
standard second-order boundary-layer asymptotic expansions for u on each of these 
sub-intervals. As u(to) is unknown, we impose the condition M as (^o) = Vo(fi;p), or 
equivalcntly u as (^o) = Vb ( — *i — £t.2+p). Here t\ and t 2 will be chosen in Lcmma [2T4l to 
match the two asymptotic expansions at x — to, while varying the parameter p will be 
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used in the construction of sub- and super- solutions. We partly follow the asymptotic 
analyses of Section 2.3.2] and [7J Section 3], where the location of the interior layer 
was outlined in the form of an asymptotic expansion t = to + st\ + e 2 t 2 + ■ ■ ■ ■ Our 
asymptotic analysis differs from these earlier works in that we expand about the 
point to instead of about the point t (which is a priori unknown); this is useful in the 
subsequent numerical analysis because our layer-adapted mesh will be centred on the 
known point to. 

As Mas(io) = Vo{0;p), the resulting asymptotic expansion Uas(x) = Un S (x;p) will 
also involve the parameter p as follows: 

(2.5) u as {x;p) := uo(x) + e 2 u 2 (x) + v {£;p) + evi(£;p) + e 2 v 2 (£;p). 
Here the smooth component uq + e 2 u 2 is defined by 

(2.6) u (x) := | M ^ xe \ t J ]} , u 2 (x):=u /b u (x,u ), 

so that F[u + e 2 u 2 ] = 0{e 4 ) for x G [0, 1] \ {t }. 

To describe the layer component vq + ev\ + e 2 v 2 , we make use of the auxiliary 
function 

(2.7) B(x,s) :=b(x,u (x) + s), 

which is clearly discontinuous at x = to. By (A2), we have -§^B(x, 0) = for all 
x to, which implies that 

(2.8) \£Lb(x,s)\ < C\s\ for x G [0, 1] \ {t }, seR, m = 0,1,2. 
Furthermore, we use the notation 



(2.9) to = io{x) 



to, x G [0, t ) 
t+ xe(t ,l] 



thus, e.g., u (i ) = (fi(to) and B(i ,s) — b(t ,(pi(t ) + s) for x < t , while Uo(*o) = 
ip 2 {t ) and B(i ,s) = b(t , <p 2 {t ) + s) for x > t . 

The zero-order boundary-layer function i>o(£) = Vo(£;p) is defined by 

(2.10) -■^v + B(io,v )=0, v (0 ± )=V (0;p)-u (t t ), v (±oo) = 0. 
Comparing this with (|2.1ap . we see that 

(2.11) v (ti;p) =Vo(Z;p)-u (i ). 

Higher-order boundary-layer components ui(£) = fi(£;p) and v 2 (£,) = V 2 (£;p) 
and defined by 

(2.12) [—g* + B s {io,vo)]v 1 = -ZB x (to,vo), v x (0) = «i(±oo) = 0, 

[-^r+S s (^ ) Wb)]« 2 = ^(0. waCO*) = -ua(4), « 2 (±oo) = 0, 
( 2 - 13 ) MO ■= -^Bxxiio^o) ~ €viB xa (t ,v ) - ^-B ss (i ,vo) 

-u 2 {i )[B s (i ,vo) - B s (i o ,0)}. 
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The functions Vi and v 2 depend on p since they are dchned using v (£]p). 

Note that v and v 2 have a discontinuity at £ = 0, but u +Vo = [uo[x)— Uo(^o)]+K) 
and «2 + i>2 are continuous at a; = to. Thus u as (a;;p) is continuous for x G [0, 1]. 

Given any suitable function v(x), introduce the functional 

*[«(•)] -^U 4 --^L =4! 

which will be used to match our asymptotic expansion at x — to ; see Lemma 12.41 

To establish the existence and properties of V\ and v 2 , note that (|2 . 12[) and (|2 . 13[) 
are particular cases of a general problem 

(2.14) [-^r+B.(t ,«b)]i/ = V(0 fo^£K\{0}, 1/(0*) = 4, i/(±oo) = 0, 

for which we have the following result. 

LEMMA 2.2. Let |^>(£)| < C(l + |£| fc )x(£)- Then there exists a solution v of 
problem (2~TJ\ ), which satisfies \u{£)\ < C(l + |£| fc+1 )x(£) and 

(2.15) $M =^(0-)- i/(0+) = 5 J 5y (- r^)x(O^+[^-^o + ]x'(0))- 

^ — oo 

Furthermore, if > and > 0, t/ien z/(£) > /or all £. 

Proof. The desired assertions follow from the explicit solution formula 

"(0=X(0 /V 2 W X(t)mdtd v + ^ X (0 for£<0, 
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K0=x(0 /_ X~ 2 (??) / x(*)^W*d»?+^yX(0 forC>0, 



which is obtained by variation of parameters noting that, by (|2.1a[) . (|2.2[) . the function 
X satisfies — jgzX + X-^sC^O) u o) = 0; see [SI Lemma 2.2]. Now, a calculation yields 
(|2.15l) and the other assertions. 

We shall show, e.g., that \v\ < C(l + |£| fc+1 )x for £ < 0. As it follows from 
(|2~4"|) . ([2~TT|) that C'x < v < C"x, then we have x*dt < ^7V dv = ^d^o) and 
therefore \tp\xdt < C(l + \t\ k )d(vQ). Thus, integrating by parts k times, one gets 
| J^xOM*) dt\ < C X 2 (l + \v\ k )- The desired bound follows. □ 

We shall now apply Lemma 12.21 to establish properties of Vo, v\ and v 2 . 

Lemma 2.3. For any constant t\ and t 2 in i2.1c\) . there exist solutions vo, i>i 
and v 2 of problems i2.10\) . h2.12\) and \2.13]) . respectively. The function vq satisfies 

(2.16) (sgnO-«D(O>0 and \v (0\ < C"x(0 for £eR\{0}. 

Furthermore, assuming that \t\ \ + \t 2 \ < C and \p\ < p* , for any arbitrarily small but 
fixed A G (0,7), there is a constant C\ such that 

(2.17) <C Ae -(7-A)|f| for £ eM\{0}, j = 0,l,2, fc = 0,...,6. 



Proof. The existence and properties ()2.16j) of the function vo as well the bound 
([2171) for i = 0, fe = 0, 1 follow from the observation (|2~TT|) combined with (f2~9]) 
and Lemma [XT] Next, the existence of v\ and v 2 and the bound (|2.17|1 for j = 1,2, 
k = are obtained by applying Lemma \'2.2\ to problems (|2.12[) and (|2.13[) . in which 
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the right-hand sides are estimated using (|2.8[) with m = 1,2 and also |i>o| < C"%. 
Similarly, the higher-order derivatives of i>o, t>i and u 2 all satisfy problems of type 
(12 . 14)) with various data, so the bound (|2.17|) for them is obtained by again applying 
Lemma [2^21 □ 

The main result of this section is as follows. 

Lemma 2.4. For the asymptotic expansion w as (x;p) from \2. 5\) we have 
(2.18a) Fu as (x; P) = 0(e 3 ) for x E (0, 1) \ {t Q }. 

Furthermore, there exist values of t\ and ti in i2.1c]) . independent of e and p, and 
positive constants C\, C% and e* = s*(p*) such that for all e < e* and < \p\ < p* 
we have 

(2.18b) (sgnp) • $KsO;p)] > C x e\p\ - C 2 e 3 . 



Proof. The relation (|2.18a[) is a standard outcome of the method of asymptotic 
expansions that was applied to generate the terms in (|2.5[) . 

To establish (|2.18bj) . note that (j2l))) implies $[u ] = e[yi(*o) - ¥>2 (*())]• As wc 
also have $[u ] = $[Vb] = and $[e 2 w 2 ] = 0(e 3 ), then 

(2.19) $K S ] = e[ip[(t ) - ip' 2 (t )} + e$M + e 2 $M + 0(e 3 ). 

By applying (|2. 15|) to problem (|2.12[) . we get 



J — oo 

Note that B x (x, vq) = b x (x, uo(x) + vq) + u' (x) b u (x, uq(x) + vq) for x ^ to- Combining 
this with u (*o) + «o = Vb yields S x (£ , v ) = b x {t Q , V ) + u' (i ) b u (t , V ) and thus 

n OO pO pOO 

^} = ^m Zb x (t ,V ) X d$+?4& Zb u (t ,Vo)xdZ+ ! '$$- $b u (t ,V ) X dt;. 



J — OO J —oo o u 

For the second term on the right-hand side, we note that b u (to, Vo)x = w[b(to, Vo(0)] 
so an integration by parts yields 

o M /.O 

$b u (t ,Vo) X dl; = - b(to,V )dt = - X '(0 d£ = -*(0). 



Combining this with a similar estimate for the third term, we arrive at 

/oo 
eMto,V )xd£-bi(*o)-¥>3(*o)]. 
-oo 

For the first term here, recall (|2.1c|) and therefore switch to the variable £ = £ — Fi + _p 
so that the resulting integral involves the functions Vo and X, which are independent 
of i\ and p, rather than Vo(£,]p) and x(£iP) : 

/•OO />00 

(2.21) / Zb x (t ,V ) X dZ= (i+t 1 -p)b x (t ,V (i))x(t)di=Cj I -(t 1 -p)C I . 
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Here 

/oo 
ib x (t ,v (i))x(i)di 
-oo 

is a fixed constant, independent of i\ and p, and 

/oo r<P2{ta) 
M*o> Vo(|)) X(l) « = - / M*o> «) dv > 

is a positive constant that appears in Assumption (A5). We now choose t\ := Cji/Ci 
so that fi = Cn/Ci + eti and therefore C77 — (tx — p)Cj = (p — et^)Ci. Combining 
this with (|2~2T)j) and (|2~2"T1) yields 

*M = Hoj(P - e*2)Cj - K(t ) - ^(io)]- 
Substituting this result in (|2.19[) . we arrive at 

= e^oyb - rf2)Cj + e 2 $b 2 ] + 0(s 3 ). 

Now, by applying (|3T5|> to problem (gUj) , we get $[w 2 ] = ^gy [O727 + 0(p + e\t 2 \)], 
where Cm equals the expression in the parentheses of formula (|2.15[) evaluated using 
the data of problem (|2.13[) in the case of p = and ti = ti; thus Cji is independent 
of p and e. Now choosing ti := Cm/Ci yields 

(2.22) $K S ] = ^ojeC/p + 0{e 2 p) + 0(e 3 ). 

Note that there exists C such that x(0) = x(— p + 1\ + £^2) satisfies > C" for all 
e < 1 and < p* . Thus choose Ci := \CiC so that eCi > 2eC x . Finally, by 
choosing e* sufficiently small and C2 sufficiently large so that the O terms in (|2.22|) 
satisfy \0{e 2 p)\ < eC x \p\ and |0(e 3 )| < C* 2 e 3 , we establish (f2~T8bj) . □ 

Note that the numerical solution of problem (jTTTJ) presents substantial difficulties 
and instabilities [6]. In that paper we describe a special numerical treatment for 
particularly small values of e that is based on the following result: 

Lemma 2.5. Let r = elnN for some C T > 2 and N > 2. Then the asymptotic 
expansion u as (a;;0) of K2. 5\) can be written as 



(x;0) = U{x,e) + 0(e\nN + N~ 2 ), U(x,e) 



V (^;0), \x-t \<r, 

Uq(x), \X - t 1 > T. 



Proof. Note that (|2.5[) immediately implies that u as (a;; 0) = uo(a;)+t>o(£; 0)+O(e). 

(i) Let \x - to I < r. By (EH), u as (a;;0) = F o (£;0) + [u (.t) - u (t )] + 0(e). Here, 
by virtue of (|2.6p and (|2.9p . one has |u (.t) — ito(*o)| < C\x — to| < Or. Consequently 
u as (x; 0) = Vb(£; 0) + 0(r + e), which yields the desired result. 

(ii) Now let \x - 1 | > r, i.e., |£| > r/e. Note that ([2~T6| combined with (|23]l yields 
I u | < 0"0 A e-^- A )^l. Choosing A sufficiently small so that C r (l - A/7) > 2, one 
gets e~^~ x ^ T / e < N~ 2 . Consequently u as (x; 0) = uq(x) + 0(N~ 2 +e) and the desired 
result follows. □ 
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3. Perturbed asymptotic expansion, sub- and super-solutions. For the 

numerical analysis that appears in [B] we now perturb the asymptotic expansion 
u as (x;p) of (|2.5[) as follows: set 

(3.1) (3(x) = P(x;p,p',k) :=u as {x-p)+p'[v^-p) + C ]+h 2 z{£;p). 

Clearly /3 is a small perturbation of u as when the parameters p' and h are small. 
In this definition, the parameter h is related to the mesh used in [5] as the compo- 
nent h z{^;p) is added to compensate for the principal part of the truncation error 
produced when the finite difference operators of [5] are applied to u as {x,t). The 
component + Co] is added to ensure that (sgnp') • F(u as + Co]) > 0. 

The functions «*(£) = f*(£;p) and z(£) = z(£;p) used in (|3.ip are defined by 

(3.2) [—^+B 8 (i ,v )]v* = \v \, for£eR\{0}, w»(0) = u*(±oo) = 0, 

(3.3) l-J^ + B s (io,v )}z=±-^V for£eR, z(0) = z(±oo) = 0. 

The functions v* and z depend onp since they are defined using vo(£;p) and Vo(£;p). 

Lemma 3.1. Assume that \p\ < p* for some positive constant p* . Then there exist 
solutions and z of problems i3. 2\) and \3. 3\) respectively, and for any arbitrarily 
small but fixed A 6 (0,7), there is a constant C\ such that 

(3.4) w.>0, \^\ + \^z\ <C A e-<^ A )l«l for £ e R \ {0}, k = 0, . . . , 4. 

Furthermore, there exist positive constants C\, Gi, C3 and e* = s*(p*) such that for 
all e < e* and < \p\ < p* we have 

(3.5) (sgnp) • p,j/, %)] > de\p\ - C 2 e 3 - C 3 \p'\. 



Proof. The existence and properties (|3.4p of v* and z are obtained by applying 
Lemma [2~2l to problems (|3.2f) and (|3.3[) . respectively. Furthermore, the relation (|2.15l) 
from this lemma yields the values of $[?;*] and $[z]. For $[«*], using the sign property 
of vq from (|2.16j) . we get 

/oo />0 poo 

Klx^=-^j(/ v oX d£- v oX dt) 
-00 J — 00 Jo 

= -Ho) (N(0~)] 2 + M0+)] 2 ) > -C 3 for all \p\ < p* . 

Here, in view of (|2.2[) , we used vox = 5(^0)' ano - ^(i 00 ) — 0, and it was understood 
that x(0) = x(0;p) and w (0 ± ) = w o (0 ± ;p). Next, for $[z], noting that ^V = X '" 
and x' = K)" = & (io,^o), we get 

/oo /.oo 
X"'X = / ^ ^ = S^OJ = 0. 

-oo J —00 00 

Thus $[/3] = $[u as ] and combining this with (f3T6|) and (|2.18bjl yields (|3~5j) . 

□ 

Lemma 3.2. There exist positive constants Co, C4, p'* and e* such that for all 
x e (0, 1) \ {t }, £ < £*, \p\< V*> < \p'\ < P 1 *, the function (3 of [3~1)) satisfies 

(3.7) (W) • [F0 - ~££;V ] > iCob'| 7 2 - C 4 (£ 3 + eh 2 + h^). 
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Proof. We partly imitate the analysis of [5J Lemma 3.2]. As, by (|2.18a[) . we have 
Fuas — 0(e 3 ), it suffices to estimate F\ u , where we use the notation F\^ := Fw — Fv 
for any two functions v and w. Noting that for it as of (|2.5[) we have u as = UQ+Vo+O(e), 
which, by (|2.7p . implies b u (x, u as ) = i? s (a;, Wo) + 0(e), we obtain 

^C +P '" + " 2Z = ~P'w v * ~ h % z + + ^)[^ «o) + °( £ +P' + ^ 

Next, using (|3 . 2() and p. 31) for -^rzv* and and then noting that for x ^ to wc 
have |_B s (x,wo) — -B s (£o,t>o)| — C|x — to|| u o < Ce, yields 

(3.8) F\ u u ^ p ' v ' +h2z = p>\v \ + S|U + \p> + P]0{e+p> + h 2 ). 

Similarly, as b u (x,u as + p'v* + h 2 z) = B s (x,vq) + 0(e + p' + h 2 ) and B s {x,vq) = 
B s (x,0) — \(x)\vq\, where \X(x)\ < C$ for some C5, we get 

< S+P ^ Z = F\::t p X + S: +Cop ' = <vp.m) - a^i + 0(e +P ' + ^]. 

Combining this with p.8[) and Fu as = 0(e 3 ), yields 

^ ~ TlW Vo = C oP' B *( x > °) + P'M(1 - CbA(x)) + 0(ej/ + p' 2 ) + 0(e 3 + eh 2 + h A ). 

Note that here B s (x,0) = b u (x,uo(x)) > 7 2 for x 7^ £cb by Assumption (A3), so 
C o \p'\B s (x,0) > C \p'\~f 2 . Now, choosing C = C7 1 so that (1 - C Q X(x)) > 0, and 
also p'* and e* sufficiently small so that \0(ep' +p l2 )\ < ^Co\p'\j 2 , we get the desired 
assertion (|3.7[) for some constant C4. □ 

Lemma 3.3. Let p >0,p' = C'ep for some positive constant C , and h 2 < Ce^ 
for some fixed (i 6 (0, 1]. Then there exists e* = e*(C", p) such that for the function 
(3 from h3.1\) we have 

(3.9) /3(x;-p,-p',h)<0(x;p,p',h) for x e [0, 1], e < e% \p\ < p*. 

Furthermore, for any arbitrarily small but fixed A £ (0, 7), there is a constant C\ such 
that u as from i2. 5\) satisfies 



(3.10) \(3(x- ±p, ±p', h) - u as (x; 0)| < C x (\p\ + h 2 )e~^ x ^ + Ce\p\. 



Proof. Fix h, and consider f3(x;p) := (3(x;p, C'ep, h). As /3 is continuous on [0, 1], 
to establish (|3.9[) . it suffices to show that > for x 7^ to- First, we note that 

|;«o = X, |^|<C(1 + ^') X fori = 1,2, |^«.| + |^|<C(l + |f|) X . 

The relation for vo here follows from (|2.1ip and (|2.2[) . The other relations are obtained 
by differentiating problems (|2.12[) . (|2.13p . p.2[) and p.3[) with respect to p, which yields 
four problems of the type (|2.14[) for J^t>i,2) an d J^z, respectively. By Lemma [2^21 
applied to this problems, the above estimates follow. Now a calculation, using p. II) . 
(|2.5|) and u* > 0, shows that for some constant C" wc have 

!^>x + c' e c -cv i (i + £ 4 )x. 



10 



By (|2.3|) . there exists a sufficiently large constant C such that for |£| > £* := C\ lne| 
we have (1 + £ 4 )x < C c ?,° e 1 ^, which implies > for £| > £*. Otherwise, if 
|£| <C'|lne|,thenC"e AI (l+f 4 ) < GV(1 + | lne| 4 ) < 1 , provided that e* is sufficiently 
small; thus we again get ■§-$ > 0. Thus we have proved (|3"1?]) . 

Similarly, one has |£/8| < C(l + £ 4 )x + CoeC". Combining this with u as (x; 0) = 

/3(x;0,0,h) — h 2 z = /3(x;0) — h 2 z and exponential-decay estimates from (|2.3[) and 
()3.4|) . we get the remaining desired estimate (|3.10[) . □ 

In this section we established the properties (|3.5I) , (|3.7|) and (|3.9[) of the function 
j3 = /3(x;p,p',ti) that are used in [6] to construct discrete sub- and super-solutions. 
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